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We study constrained Hamiltonian systems by utilizing general forms of time discretization. We 
show that for explicit discretizations, the requirement of preserving the canonical Poisson bracket 
under discrete evolution imposes strong conditions on both allowable discretizations and Hamil- 
tonians. These conditions permit time discretizations for a limited class of Hamiltonians, which 
does not include homogeneous cosmological models. We also present two general classes of implicit 
discretizations which preserve Poisson brackets for any Hamiltonian. Both types of discretizations 
generically do not preserve first class constraint algebras. Using this observation, we show that time 
discretization provides a complicated time gauge fixing for quantum gravity models, which may be 
compared with the alternative procedure of gauge fixing before discretization. 

I. INTRODUCTION 

In all approaches to quantum gravity the idea of a discrete structure of spacetime at the Planck scale is fundamental. 
In lattice approaches, such as Regge calculus, dynamical triangulations, and other path integral methods, discreteness 
in space and time is built in right from the start 0. In canonical methods, such as loop quantum gravity, discreteness 
of space arises as a rather natural consequence of Dirac quantization 0, OJ 0, HJ , and it is possible that a form of 
discrete time evolution may emerge from a regulated action of the Hamiltonian constraint. This is already manifested 
to some extent in spin foam models 6], which may be viewed as " covariantizations" of loop quantum gravity models. 

In attempts to construct a path integral for quantum gravity, a standard starting point is a classical non-Hamiltonian 
discretization using triangulations. It is interesting to ask whether a phase space path integral for quantum gravity 
can be defined starting from a classical discretization that is Hamiltonian. If possible, this would provide a reg- 
ularization that is more closely tied to the classical theory than triangulations or spin foam methods. There are 
two main approaches to the problem: (i) fix all coordinate gauge conditions classically and discretize the reduced 
Hamiltonian system, or (ii) attempt a discretization of the constrained system and study its consequences for the 
Hamiltonian structure. The primary requirement for canonical quantization of a discrete system is to ensure that the 
fundamental Poisson bracket is preserved under discrete classical evolution. In the second approach, the consequences 
of discretization on the constraint algebra also require scrutiny. 

These questions have been studied to some extent in the literature. Discrete evolution schemes for unconstrained 
Hamiltonian dynamics appear in Q, with further refinements in where schemes preserving the Poisson bracket 
are presented. More recently, the applicability of discrete evolution schemes to constrained Hamiltonian theories have 
been discussed Q , with a view to applicability to quantum gravity. 

From elementary numerical methods, it is known that there are implicit and explicit approaches to discretization 
of differential equations. Although cither approach may be useful for studying classical evolution, explicit methods 
are more natural for quantization. This is because discrete Hamiltonians that are functions of canonical variables at 
more than one time instant are effectively non-local in time, so it is unclear whether canonical quantization makes 
sense. It is therefore necessary to study classical discretizations for which the discrete Hamiltonian does not depend 
on more than one time point. This may restrict somewhat the potential applicability of implicit schemes. 

For the purpose of developing a fundamental lattice approach to quantum gravity starting from a discretization of 
Hamilton's equations, it is necessary to ensure that canonical structures are not affected by the procedures used. This 
is because classical preservation of Poisson brackets is a requirement for a clean approach to canonical quantization. 
The main purpose of this paper is to investigate this aspect. We study certain issues associated with discrete time 
evolution. We do not address space discretization, since our main concern is the extent to which time discretization 
can maintain Hamiltonian structures for quantization. Specifically we address the following questions: (i) What are 
the general classes of discretization? (ii) What are the conditions required to maintain the fundamental Poisson 
brackets under evolution for a given Hamiltonian system? (iii) Are there restrictions on allowed Hamiltonians? (iv) 
How does the approach work in simple systems, and (v) How viable is the approach for quantum gravity? 

Our starting points are general forms of explicit and implicit discretizations containing arbitrary functions of phase 
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space variables. In Section II we derive the necessary conditions on these functions to preserve the fundamental Poisson 
bracket for both types of schemes. We also study the consequences of discretization for constraint evolution for the 
parametrized particle. In Section III we apply the approaches to the parametrized scalar field, FRW cosmological 
models, and BF theory. In the final section we dicsuss some consequences of our results for applicability to quantum 
gravity. 



II. DISCRETE HAMILTONIAN EVOLUTION 



As for differential equations, there are two general approaches that may be utilised to obtain a discretization of 
Hamilton's equations that preserve the canonical Poisson bracket. The explicit schemes, as the name suggests, give 
explicit functions of the form 

Pn+i = f(q n ,p n ), q n +i = g(q n ,p n ) (l) 

for evolution. This is in contrast to implicit schemes, where the general equations are of the form 

F(p n+1 ,q n+1 ,p n ,q n ) = 0, G(p n+1 ,q n+1 ,p n ,q n ) = 0. (2) 

We consider both types of discretization, and outline the conditions under which the schemes preserve Poisson brackets, 
and thus have a utility for quantization at the discrete level. 



A. Explicit schemes 



Consider a one particle system with canonical phase space variables (q,p), and arbitrary Hamiltonian H(q,p). A 
first attempt at time discretization of Hamilton's equations with step e might be the simplest finite differencing scheme 



q n +i = q n + t 



Pn+l = Pn~ e 



dH(q n ,p n ) 

dp n 
dH(q n ,p n ) 



dq„ 



(3) 



However there is an immediate problem with this: it does not preserve the Poisson bracket under evolution for even 
the standard Hamiltonian H(q,p) = p 2 /2m + V(q); the order e 2 term does not vanish unless V(q) = 0. To rectify 
this it is necessary to introduce higher order terms in e in the defining equations of the discretization. This might 
be done using explicit or implicit finite differencing schemes. We choose a general explicit approach and study its 
consequences. The reasons are simplicity and potential utility for canonical quantization. 
The defining equations for the discrete evolution we study are the following: 



q n +i = q n + £ 



Pn+l 



Pn 



dH(q n ,p n ) 

dp n 

dH(q n ,p n ) 
dq n 



+ e 2 a(q n ,p n ) 
+ e 2 P{qn,Pn)- 



(4) 
(5) 



a, (3 are at this stage arbitrary functions. The Poisson bracket of the evolved variables is 



{q n +i,Pn+i} = 1 - e 



oq n J 
{a,p n } + {<?„,/?} 



dp n J 

J ' dH dH \ 

\ dp n ' dq n J 



{ dp n ' ^ 



' dq n J 



+e 4 {a, (3} 



(6) 



This equation displays the conditions required for preserving the canonical Poisson bracket. The first order term in e 
vanishes identically because the Hamiltonian is a function of the discrete phase space variables at a single time point. 
(See below for discussion of a discretization where this is not the case.) The remaining task is to see what functions 
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H, a, and (3 lead to vanishing of the three higher order terms in e. It is obviously better to fix H by looking at 
systems of interest, and then attempt solutions for a and (3. 

In order to obtain a general procedure for solving these conditions, it is useful to consider first the Hamiltonian 
independent e 4 term. Its solutions can be classified into one of the following four cases: (i) a = 0, (ii) (3 — 0, (iii) 
a(q n ), (3{q n ) or a(p n ), f3(p n ), and (iv) a = (3 ^ 0. For a given Hamiltonian, the other conditions can be written out 
explicitly. 

The restrictions that arise may be seen by considering the standard Hamiltonian 

H(q,p) = ^ + V(q) (7) 
The conditions have a unique solution: Case (ii) gives 

«.) - ^ (8) 
m oq 

The remaining cases (i),(iii) and (iv) do not give a solution. This result is surprising given that there are numerous 
possibilities for discretizing classical differential equations. The lesson is that preserving canonical Poisson brackets 
imposes rather strong conditions on allowable discretizations. 

As a further illustration consider the parametrized particle. This case is also useful for seeing what happens to the 
constraints under discrete evolution. The canonical variables are (i,p 4 ) and (q,p), which depend on a time parameter 
A. The Hamiltonian constraint is 

H := p l + |i + V(q) = (9) 
Following (JSJ, discrete evolution with lapse function N is defined by 



<7n+l 


= q n 


+ Ne 


Pn i 2 / \ 
he a(q n ,p n ) 

m 


Pn+1 


= Pn 


-Ne 


dV{q n ) t2a , ^ 
„ +e P{q n ,Pn) 
dq„ 


Pn+1 


= Pn 






tn+1 


— t'n 


f Ne. 





(10) 

From these equations it is clear that preservation of the Poisson bracket goes through in the same way as for the 
unparametrized case, with f3 = and a given by JHJ up to a factor of N. 

The second check is the evolution of constraints using Q10[l. Is H n +i = i?(t n +i)P^ + D In+liPn+l) = H n l A 
straightforward computation using the evolution equations gives 

H n +i =p t n + 1 L [p n -eN8V + e 2 f3] 2 + V(q n+1 ) 

2 

- H n + i- [N 2 {dVf +/3 2 + 2 Pn f3 + 2ma] + 0(e 3 ), (11) 

where dV = dV(q n )/dq n . From this it is clear that H n +i does not equal H n : the Taylor expansion of the potential 
contains terms that cannot be cancelled by any choice of the functions a, j3 (which are in any case effectively fixed 
by the requirement of preserving the canonical Poisson bracket). The discretization has made the constraint second 
class. 

At this stage there are two possible ways to proceed: (i) impose the Hamiltonian constraint strongly since it is not 
preserved in time, or (ii) fix the free function N such that the additonal terms on the right hand side of i|ll|) vanish. 
The latter amounts to N becoming a fixed function of n, and hence the time step eN also becomes a function of n. 
Thus, making the evolution consistent via (ii) requires a variable time step in this sense, which concomitantly means 
fixing the lagrange multiplier. 

From the theory of constrained Hamiltonian systems it is possible to see that these two choices are in fact equivalent: 
Gauge fixing leads to first class constraints becoming second class pairs with the gauge conditions. Since we started 
with a first class system (the parametrized particle), which became second class as a result of time discretization, 
we can ask what time gauge condition has been implicitly imposed by the discretization that led to the Hamiltonian 
constraint not being preserved in time. 
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To answer this question let us recall the procedure for the continuous case in the canonical theory, where time 
gauge fixing is accomplished by setting a function of the phase space variables to equal the time parameter. Consider 
a phase space function /(i,p 4 , q,p) and impose the gauge condition 

X :=A-/ = 0. (12) 

This fixes parameter time A. Preserving the gauge condition under evolution requires 

! = | + {x ,iVIO = 0. (13) 



This fixes the lapse as a function of / 



N = (14) 

{f,H} K > 



Conversely, fixing N first leads via the above equations to a connection with the canonical gauge fixing function /. 
Thus a discretization may also be determined by fixing the lapse, which may be viewed as imposing a variable time 
step. This of course, is what is routinely done in numerical relativity, with the difference that in that case it is not 
necessary to preserve the canonical Poisson bracket for quantization. 

For the discrete evolution we have defined, this procedure may be used to extract the gauge fixing condition. This 
is done by first imposing H n+ i = strongly since it is second class. This gives an equation for N. The gauge 
fixing function / is then determined by solving (|14f) . which in general is rather messy. For example, for the above 
discretization of the parametrized particle, it is not possible to solve for N to all orders in e due to the presence of all 
orders of N in H n+ \ (|llfl . (This is because the expansion of V(q n +i) in powers of Aq n = q n +i — q n contains powers 
of N.) Thus only an approximate solution for N is possible, up to a specified order in e. 

The lessons from this section are two-fold: (i) The conditions on allowable explicit time discretizations are rather 
strong if canonical structures are to be preserved, (ii) For constrained theories with time reparamctrization invariance, 
time discretization is equivalent to rather complicated gauge fixing conditions. 

B. Implicit schemes 

There are implicit schemes for time discretizations discussed in the literature. We mention two examples below 
before introducing somegeneralizations. One of these uses a leap-frog numerical scheme to construct a consistent 
Hamiltonian evolution The other method is implicit and uses a " Hamiltonian" that depends on phase space 
variables at neighbouring time steps [9|- 

The leap-frog discretization scheme for the standard Hamiltonian p 2 /2m + V(q) is Q 

p n e 2 dV{q n ) 

q n +i = q n + e „ 5 

m 2m oq n 
e f dV(q n ) dV(q n+1 ) \ 

Pn+l = Pn - ~ — g + — £ ■ (15) 

2 V dq n oq n+ i J 

This is an implicit scheme since the V(q n+ i) term in the p n +i equation is determined by first computing q n +\. 
Without this term it is, up to factors, the same as the one we derived above from the general form (J5J). It preserves 
the canonical Poisson bracket under evolution. 

The non-local Hamiltonian scheme is rather unusual. Evolution equations are derived from a discrete action, 
which may be viewed as a function of coordinates q n and velocities v n : 

S[q n ,v n ] = y^^L n (q n ,v n ) = ^ [v n (q n +i - q n ) - eH(v n ,q n )} . (16) 

n n 

This resembles a Hamiltonian form. The conjugate momentum is defined by p n +\ '■= dL n / dq n +\ = Vy. This leads to 
the following discrete equations for the canonically conjugate pair (q n ,Pn) (Eqns. (37,43-44) of Ref. yj): 

dH(q n ,p n+1 ) 

q n +i = q n + i- — k 

OPn+1 

dH(q n ,p n+1 ) 

Pn+l = Pn~ £ 5 (17) 

oq n 
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It is in the sense manifested by these equations that the discretised phase space Hamiltonian is non-local in time. 

The requirement of preserving the Poisson bracket for more general implicit schemes is at first sight more complicated 
than the procedure we followed for the explicit case. The computation of Poisson brackets requires implicit differentials 
assuming the validity of the implicit function theorem for the discretization. For a given discretization this leads to 
restrictions on the form of the Hamiltonian. It therefore appears that a general result may be difficult to obtain, and 
that one is stuck with a case by case analysis. 

There is however a general approach for finding a large class of implicit schemes that are automatically consistent. 
This involves (i) viewing discrete evolution as a canonical transformation from the variables (q n ,Pn) to the variables 
{ln+i>Pn+i)> an< i (h) requiring that the canonical transformation resemble a discretization of Hamilton's equations. 
This approach gives two general classes of implicit discretization starting from a given Hamiltonian H(q,p). 

To see how this approach works it is useful to consider each of the four classes of generating functions for canonical 
transformations. As a first example consider a generating function of the form 

I(q n +i,p n ) = -q n +iPn + e F(q n +i,p n ) + e 2 G(q n+ i,p n ) H (18) 

The standard canonical transformation rules give the equations 

qn+i = 

Pn+i = Pn-e — — e 2 1 . (19) 

dq n+ i oq n+ i 

The discretization prescription is to fix the functional form of F to be that of the given Hamiltonian H(q n ,p n ) with 
the replacement q n — > q n +i, ie set 

F = H(q n+1 ,p n ). (20) 

The result is an automatically consistent class of implicit discretizations where the function G is arbitrary. The limit 
e — ► gives the continuous time Hamilton equations. 

The second example of this type is provided by the generating function 

I{q n ,Pn+i) = -q n Pn+i + e F(q n ,p n+ x) + e 2 G(q n ,p n+1 ) H (21) 



dF 


2 dG 


OPn 


OPn 


dF 


2 dG 


dq n +i 


dq n+ i 



The transformation equations are now 



dF 2 dG 



q n +i = q n + e he 

opn+i dp n+ i 

dF 2 dG 

Pn+i = Pn-e q e 1 (22) 

dq n oq n 

These again give the continuum Hamilton equations in the e — > limit with F — H(q n ,p n +i) for arbitrary G. (A 
special case of this scheme arises from the action viewpoint discussed in |9j.) 

The remaining two types of transformations with generating functions depending on the pairs (q n iQn+i) an d 
iPniPn+i) do not give a form resembling discretized Hamilton equations, so the two cases discussed above appear to 
exhaust this approach for generating consistent implicit schemes. 

As for the explicit scheme of the last section, it is readily verified that constraints are not preserved generically 
with these two types of implicit discretization. For the parametrized particle @ for example, the evolution equations 
with the first method (|T$)l (without the e 2 terms, and with lapse N) are 

Pn 



q n+ i = q n + Ne 



p n+ i = p n - Ne 



t _ t 

Pn+l — Pn 



m 

dV(q n+1 ) 



dq n 



+i 



t n +i = t n + Ne. (23) 

Substitution of these into H n = p n + p 2 /2m + V(q n ) shows that H n +i ^ H n . 

In summary, we have seen that both the explicit and implicit schemes are quite restrictive in their respective 
settings, although with the latter applicable to any Hamiltonian. In particular, the schemes we have considered are 
two-step, involving only variables at times n and n+1. Multi-step general schemes cannot be represented as canonical 
transformations in any simple way for the implicit case, and checking of Poisson bracket preservation appears to be 
much more tedious for the explicit case. It is also apparent from the simple example of the parametrized particle that 
constraints are not preserved under discrete time evolution. 
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III. APPLICATIONS 



It is useful to see whether there are viable implementations of the explicit and implicit discretization discussed in 
the last section. In the following we consider a few examples of models with Hamiltonian constraints to see if there are 
solutions for the functions a, j3 in iJSJ that preserve the canonical Poisson bracket. To do this we follow the procedure 
described above. We also examine the time gauge fixing conditions introduced by the implicit discretizations I|19I22[1 . 



A. The parametrized scalar field 



As a first application we look at the parametrized scalar field in a curved spacetime. Its Hamiltonian is 



H 



N 



V 1 + 27/| H 2 + \ Vq q ah d a ct>d b <f> + sfq V{cj>) 



(24) 



where the canonical coordinates are and (</>, II), and q a b is the spatial background metric. The discrete evolution 

equations are 

4> n+ i = 4> n + eN -^- + e 2 a(<f) n ,n n ) 

n n+1 - n„ + e [d a (^q ab d b ^N) - V'Ny/q] + e 2 /3(0„,II n ) 

t n +i = t n + Ne 



Pn+l 



(25) 



From this one calculates {0„+i, II n +i} to try to find functions a and f3 such that {<p n +\, II„+i} = 1. It turns out that 
this is not possible except for the special case (f> = <fi(t), where it reduces to the case of the parametrized particle 10. 



B. Cosmological models 

We consider two models here, FRW coupled to a scalar field and deSitter space. These serve to illustrate the main 
problem with discretization arising from the Hamiltonian constraint of general relativity: the kinetic term contains a 
mixture of gravitational configuration and momentum variables, which makes it difficult to find a discretization that 
preserves the canonical commutation rules under evolution for the explicit scheme. 

For both models we use geometrodynamics phase space variables (q a fc, 7r ab ), and the parametrization 

q ab = a(t) e ab n ab = i P(t) e ab (26) 

where e ab is the flat Euclidean 3— metric. (The main features of the analysis are unchanged for the triad-connection 
phase space variables.) 

The Hamiltonian constraint for flat FRW minimally coupled to a scalar field is 

H(a,p) = -\ a}'V + \ a" 3 / 2 H 2 a 3 / 2 V(ct>), (27) 

where V((f>) is a potential for the scalar field 0, and n is its canonical momentum. For De Sitter space, the potential 
is replaced by the cosmological constant A. 

The discrete evolution equations obtained using JSJ are 

a„+i = a„ - eN i al/ 2 p n + e 2 a(a n ,p n ) (28) 
for both the FRW and deSitter case. The momentum evolution equation is 

+ e 2 (3(a n , Pn ) (29) 
for the FRW case, and 

+ e 2 /3(a„,p„) (30) 



Pn+l = Pn 



eN 



12 



a _1/ V 



ii' 2 iv 2 



,1/2 



Pn+l = Pn + eN 



y 2 a- 1/2 pl + \ 4 /2 A 
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for deSitter. (We ignore the evolution equations for the scalar field, since these are not relevant for the point we wish 
to make.) The main observation again is that computing {a„ +1 ,p I1+1 } shows that there is no choice for a(a n ,p n ) and 
/3(a n ,p n ) which preserves the Poisson bracket to all non- vanishing orders in e. 

Thus we conclude that this general form of discretization is not applicable to models derived from general relativity. 



C. BF theory 



This example of a fully constrained dynamics is of particular interest in the context of spin foam models. Once 
again we focus on time discretization, which is what is relevant for preserving the Poisson bracket under evolution. 



Abelian case : the continuum Hamiltonian is given by 

H = J (Xe ab F ab + [id a E a ) d 3 x 



(31) 



with the two constraints F ab — and d a E a — and Lagrange multipliers A and /i. The discrete evolution equations 
for the canonical variables A a and E a are 

A a ,n+i = A an — ed a ^n, (32) 

K+i = EZ + e2e ab d b \ n . (33) 

From these it is obvious that the Poisson brackets are trivially preserved, even without considering terms of order e 2 . 

Non- Abelian case: The only difference is that the canonical variables are now Lie-algebra valued. Thus the Hamiltonian 
is 

H = J {X l t ab F l ab + ^D a E al )d 3 x (34) 
where D a = d a + A a is the covariant derivative. As a result, the discrete evolution now becomes non-trivial: 



4,n+i = 4,n-e D af in + e 2 a(Al n ,E™) 



(35) 
(36) 



As a consequence of the additional terms compared to the Abelian case, it is straightforward to check that there is 
no solution for a and f3 that leads to a preservation of the Poisson bracket. 

For the implicit schemes the Poisson bracket is of course preserved by the general argument given in the last 
section. The constraint associated with time evolution in BF theory is a combination of the two " internal symmetry" 
constraints, as occurs for example in 2+1 gravity 10]. Thus it is expected that there is some partial breaking of these 
symmetries due to time discretization. Given this, it is useful to see what are the (partial) gauge fixing conditions 
induced by the time discretization. 

Denoting the continuum gauge conditions by f l (A,E) = and g z (A, E) = 0, the relation between these and the 
lagrange multipliers A* and (J, 1 are given by the requirement that the gauges are preserved in time. Schematically this 
is 



{f,F[X] + G\n]} = 0, 



{.g,F[A] + G[ M ]} = 0. 



(37) 



For the discretization, one can in principle extract the functions X 1 and fj, 1 by setting the time evolved constraints to 
zero, ie. G l n+1 = G l (A n+1 ,E n+ i) — and F^ +1 = F l (A n+ i, E n+1 ) — 0. One can then attempt to solve Eqn. (|3T|) for 
the gauge conditions / and g. To get an idea of what the resulting equations look like, the first task is to extract A* 
and \i l from the evolved constraints. For the first implicit scheme (|19fl for example, where H = H[A n +±, E n ], these 
equations are 



n+1 ~ u 

pi jpi 

r n+l — r n 



-8- 



SH 



5A 



n+l 



ijk / pak ° H _ A3 



SEi 



SH 



_ e 2 e ijk 



5H SH 

5Ei 8A* +1 







n SH JJJ, . A 8H 

d m +2e Ai m 



SH SH 

SEl SE k n 



(38) 
(39) 
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up to order e, where we have suppressed e and the spatial indices (a, b ■ ■ •). Inserting 

SH 



— = e^4 +1 /4 (40) 



SEi 



and 

SH 



SAL 



= 2 (dXi + ( '' i ,l;,. : A!;) (41) 



into the equations above, one can in principle solve for A and /i, and substitute in (|37f) to get the functions / and g. 
These equations give an indication of the rather complex gauge fixing induced by this discretization. 

Note that although we are dealing here with time discretization only, this may still affect the first class nature of the 
purely internal symmetry constraints: the first class property is preserved automatically under evolution if G^ +1 = 
follows from G % n = 0, without fixing any Lagrange multipliers. We have seen above that this is not the case. Therefore, 
although time discertization does not affect the internal symmetry algebra at fixed time in), it nevertheless makes 
the constraints second class with respect to evolution in time. 

From the above examples, the pattern of calculation is clear for application to quantum gravity. It is unlikely that 
explicit schemes will work due to non-conservation of the Poisson bracket. The two classes of implicit schemes can 
be pursued to some extent at least, in the connection-triad canonical variables. A first attempt might be to keep 
space continuous and discretize time for the general form of the constraints with the Barbero-Immirzi parameter. It 
is possible there are choices of this parameter that simplify the form of the constraints and gauge fixing conditions. 



IV. DISCUSSION 



We have studied general classes of explicit and implicit time discretizations of Hamilton's equations containing 
arbitrary functions. The general forms were motivated by simplicity, and relative ease in establishing the primary 
criterion that the Poisson bracket is preserved under discrete evolution. 

The explicit class of schemes we studied preserves the Poisson bracket for only a limited class of Hamiltonians. 
Therefore its use for quantum gravity is rather limited. The two classes of implicit schemes have more potential 
utility in that these apply to any Hamiltonian. However the resulting gauge conditions are quite complicated. 

Our approach appears to exhaust explicit and implicit two-step discretization schemes that preserve the Poisson 
bracket. In both cases, continuum first class constraints are not preserved in time. Therefore gauge fixing conditions 
are implicit in the discretization, and can in principle be explicitly recovered. Given this one can ask what is the 
utility of this lattice approach to gauge fixing, versus the alternative path of gauge fixing in the continuum theory 
followed by discretization (where preservation of constraints is obviously not an issue). For example, the problem of 
time is "solved" by gauge fixing in both approaches. Furthermore, at the level of quantization, the old question of 
whether different gauge fixings lead to the same quantum gravity model would reemerge in a lattice disguise. 
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